hep-ph/9512208 



CERN-TH/95-323 
Edinburgh 95/560 



UNIVERSALITY AND SCALING 
IN PERTURBATIVE QCD AT SMALL x 



Stefano Forte^")* and Richard D. Ball^^'^H 

Theory Division, CERN, 
CH-1211 Geneve 23, Switzerland. 

Department of Physics and Astronomy 
University of Edinburgh, Edinburgh EH9 3JZ, Scotland 



Abstract 

We present a pedagogical review of the universal scaling properties displayed 
by the structure function F2 at small x and large as measured at HERA. We 
first describe the derivation of the double asymptotic scaling of F2 from the 
leading-order Altarelli-Parisi equations of perturbative QCD. Universal next-to- 
leading order corrections to scaling are also derived. We explain why the universal 
scaling behaviour is spoiled when the initial distributions rise too steeply by con- 
sidering the nonsinglet distribution F2 — F2 as an explicit example. We then 
examine the stability of double scaling to the inclusion of higher order singulari- 
ties, explaining how the perturbative expansion at small x can be reorganized in 
such a way that each order is given by the sum of a convergent series of contribu- 
tions which are of arbitrarily high order in the coupling. The wave-like nature of 
perturbative evolution is then shown to persist throughout almost all the small x 
region, giving rise asymptotically to double scaling for a wide class of boundary 
conditions. 
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1. Structure Functions at Small x 

The small x limit of structure functions measured in deep-inelastic scattering is the 
frontier of perturbative QCD both from the experimental and the theoretical point of 
view. In deep-inelastic scattering experiments [0, one measures the total, fully inclusive 
cross section for scattering of a virtual photon with virtuality —Q^ over a nucleon. The 
center-of-mass energy of the collision is given by the Mandelstam invariant s = Q^^—^. 
In the large limit the cross section is parametrized by a single form factor F2{x,Q'^) 
which is determined by the underlying partonic degrees of freedom, because its moments 
with respect to x are directly related to the nucleon matrix elements of quark and gluon 
operators. The scale dependence of these moments is governed by renormalization group 
equations which summarize the dynamical content of perturbative QCD. 

The small x limit at large thus corresponds to probing the light-cone dynamics 
of the nucleon (large Q^) in the high energy limit (large s ~ Q^/x). This limit stretches 
perturbative QCD towards its nonperturbative frontier. On the one hand, one might expect 
the high-energy limit of total cross sections to be governed by unitarity, and in particular by 
the t-channel exchange of Regge trajectories , which do not admit a simple perturbative 
interpretation. On the other hand, perturbation theory itself in this limit is somewhat 
problematic, hinting to its eventual breakdown. Indeed, the anomalous dimensions which 
govern the perturbative scaling behaviour grow without bounds, implying that Bjorken 
scaling is shifted to larger and larger values of Q^. This is a manifestation of the fact 
that there is now another large scale in the theory besides Q^, namely s itself: hence the 
renormalization group will have to be adapted in order to sum up this scale too. It is not 
a priori obvious in which kinematic region this might be possible, if at all. 

Experimentally, accessing this region requires the very high center-of-mass energies 
that have only been attained very recently at the electron- proton collider HERA @-|^. 
When the first data on F2 at large and small x were first presented they have provoked 
a considerable amount of surprise: they not only displayed the expected large violations 
of Bjorken scaling, growing larger as x decreases, rather, they also deviated from the 
Regge behaviour which is well tested by high-energy elastic scattering data, by displaying 
a marked rise as x decreases at fixed Q^, whereas Regge theory would have a fiat or almost 
fiat behaviour. Yet, such a violation of Regge behaviour was predicted more than twenty 
years ago as a direct consequence of the leading-order renormalization group equations of 
perturbative QCD 0. 

This non-Regge rise takes the form of a simple universal scaling law satisfied by 
F2(x,(5^) at large and small x: the structure function depends only on a variable 
a(x, Q^). Furthermore, this dependence is universal, unlike Bjorken scaling, where asymp- 
totically structure functions depend only on x, but in a non- universal, uncalculable way: 
hence it actually corresponds to a double scaling law, as the universal dependence may be 
scaled out. The way this double asymptotic scaling behaviour has arisen out of the HERA 
data is shown in fig. 1, which displays F2 as a function of the scaling variable a"(x, Q^) (for 
all accessible values of x and Q^), for each successive published set of data. 

It is apparent that double asymptotic scaling is the foremost feature of F2 at small 
X and large Q^, at least in the region presently explored by the HERA experiments. 
Understanding perturbative QCD in this region thus means understanding the physics of 
double asymptotic scaling. This entails understanding why the simple scaling prediction. 
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which, after all, follows from a leading-order renormalization group analysis, survives the 
problems of perturbative instability and need for the inclusion of other large scales alluded 
to above. 

In these lectures we will briefly summarize the current status of the current theoretical 
and phenomenological understanding of double scaling. In sect. 2 we will derive the dou- 
ble scaling prediction from leading-order perturbative QCD and discuss the physics behind 
this behaviour. We will then see how double scaling may be spoiled if parton distributions 
at low scale are too steep, taking as an example the case of nonsinglet structure functions, 
where this actually happens. We will further explain how double scaling is modified (but 
its universality preserved) by the inclusion of next-to-leading order corrections and show 
that present-day data are perfectly described by such a next-to-leading order analysis. 
In sect. 3 we will then discuss the behaviour of the perturbative expansion of anomalous 
dimensions when the effects of the other large scale which is present in the problem are 
included to all orders in the coupling. We will show that the perturbative expansion may 
be reorganized so as to sum up these effects, and that double scaling emerges then as the 
generic universal asymptotic behaviour. We will conclude by recalling some recent phe- 
nomenological applications of this formalism to precision tests of QCD, and summarizing 
the most promising future theoretical and phenomenological developments. 



2. Double Asymptotic Scaling at Leading and Next-to-Leading Order 

The perturbative evolution of the structure function F2(x, Q^) is determined by first 
decomposing it into parton distributions: 



x-^F^ix, Q2) = J2 e^C, ® {q, + q,) +Cg®g, (2.1) 



where n/ is the number of active fiavors, is the electric charge of the quark distribution 
qi{x] Q^), ® denotes the convolution with respect to i.e. [/ ® = f 9{y): 
and the coefficient functions at leading order are simply 

a{x,Q^)=6{l-xy, Cg{x,Q^) = 0, (2.2) 

while at higher orders they depend on the specific choice of renormalization prescription 
(factorization scheme): in most of the subsequent treatment we will choose a scheme 
(parton scheme) in which eq. ( |2.2| ) remains true to all orders. The evolution of parton 
distributions is in turn determined by the Altarelli-Parisi equations 

2tX \Fqg FqqJ \qSy 



d 




di 
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where t = ln((5^/A^) and the singlet and nonsinglet quark distributions are respectively 
given by 



(2.4) 



with (e^) = ^ Yl^=i ^i- terms of these distributions the decomposition (|2.1|) becomes, 



in a parton scheme ( |2.2| ), simply 

F2(x,Q2) = {e^)x{qs{x,t)+q^s{x,t)). (2.5) 

The splitting functions depend on t through the strong coupling as(t), and at n-th per- 
turbative order are given by 



27r 



P(x,t) = X:i^«(^)VoK+^). (2.6) 

i=i V ^ / 



The Altarelli-Parisi equations eq. (|2.3| ) are simply the inverse Mellin transforms of 
the renormalization group equations satisfied by the nucleon matrix elements of quark and 
gluon operators, namely 

d ( g{N,t) \ _ a,{t) (^gg{N,as{t)) 'yg,{N,as{t))\ ( g{N,t) 



|gNs(iV,t) = ^7,T(^'«^W)5Ns(iV,t). 

The operator matrix elements are related to the corresponding parton distributions by 
Mellin transform, i.e. the matrix element of the (leading twist) spin + 1 operator is the 
A^-th moment of the corresponding parton distribution 

p{N,t) = M[p{x,t)]= [ dxx^p{x,t), (2.8) 

Jo 

where p{x, t) is any linear combination of parton densities. Likewise, anomalous dimensions 
are found by taking moments of the splitting functions 

'y{N,asit)) = M[Pix,t)]= [ dxx^P{x,t). (2.9) 

Jo 



Notice that only operators for odd values of A^ > in eq. (|2.8|),(pT9|) actually exist, because 
even (odd) spin operators are charge-conjugation even (odd), while F2 and the parton dis- 
tributions which contribute to it are charge- conjugation even; the anomalous dimensions 
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for all other (complex, in general) values of can only be defined by analytic continua- 
tion. This continuation is provided by the Altarelli-Parisi formalism, where the primary 
quantities, namely the parton densities and splitting functions, can be used to define the 
values of matrix elements and anomalous dimensions for all through (|2.8| ) and ( |2.9|) , 
provided only that they are known for all values of < x < 1. Of course, these are 
also precisely the quantities which can be extracted directly from the structure functions 
measured experimentally, albeit in practice only over a limited range of x. 

The Altarelli-Parisi equations determine the parton distributions at (x', t') in terms 
of their values for all x > x' and t < t' , hence they actually describe evolution in the whole 
{x,t) plane, even though the renormalization group equations eq. ( |2.7| ) for each value of 
N only specify an evolution law with respect to t. The evolution with respect to x, while 
causal, is however nonlocal. The basic idea behind double scaling is the realization that at 
small X the Altarelli-Parisi equations actually reduce to local evolution equations in both 
variables, x and t, which can then be treated symmetrically. 

This can be shown by constructing a systematic approximation |^ to the evolution 
equations eq. ( |2.3| ) at small x and large Q^. As x gets smaller, one would expect the 
behaviour of parton distributions to be dominated by that of their Mellin transforms 
for small N, and thus by the small behaviour of the anomalous dimensions. Since the 
anomalous dimensions are singular at small N one would specifically expect their rightmost 
singularity to provide the dominant behaviour. 

A simple way of showing that this is indeed the case is to solve the Altarelli- 
Parisi equations by Mellin transformation, i.e. solve the renormalization group equations 
eq. ( p.7|) for the eigenvectors of the anomalous dimension matrix, which are linear combi- 
nations p{x, t) of the parton densities. Using the leading order (LO) form of the anomalous 
dimensions (which is t independent) and of as = ^ with (3o = 11 — ^n/, this procedure 
gives generically 

p{N,t)=p{N,to)exp ^ 



Po 



(2.10) 



where 



The X space solution is then the inverse Mellin transform 

/ioo 
dN ex.p{^N)p{N,t), (2.12) 
-ioo 



where 

e^ln(^) (2.13) 

and the integration runs over a contour located to the right of all singularities of 7(A^) and 
p{N, to). Assume for the moment that any singularities of the initial condition p{N, to) are 
to the left of those of 7(A^) (which are always poles on the real axis at non-positive integer 
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values of A^). Then at small x, i.e. as ^ grows the integral may be evaluated by the saddle 
point method: the saddle point condition is of the form 



6 + 



2 c^7, 



(1) 



(3q dN 



^C = 0, 



^2.14) 



where "f^l^g is the leading singularity of the one loop anomalous dimension. Higher order 
terms in the expansion around the singularity then give subleading corrections both to the 
location of the saddle point and to the integral over it. If the initial condition QnsI-^j^o) 
has a singularity (typically a branch point) to the right of that of the anomalous dimension 
this will also contribute to the integral and indeed may dominate the contribution from the 
saddle point as ^ — oo. However the saddle point always gives the dominant contribution 
at large i.e. large C, and in any case the dominant contribution to the evolution is 

always given at large ^ by the leading singularity j^l^g- 

We may thus obtain the leading small x behaviour of the Altarelli-Parisi equations 
by expanding the matrix of anomalous dimensions around its rightmost singularity, deter- 
mining the corresponding splitting functions by inverse Mellin transformation, and then 
solving the resulting simplified Altarelli-Parisi equations. The singularity is located at 
iV = in the singlet case, and at = — 1 in the nonsinglet. This means that all other 
things being equal q]<(s{x,t) will display the same qualitative behaviour as xqs{x,t), and 
xg{x,t), i.e. at small x the nonsinglet will be down by a power of x compared to the 
singlet. 

Expanding the singlet anomalous dimension about its leading singularity, at one loop 
we find 



1 



2Ca 




2Cf 




+ 
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+ OiN), 



while doing the same thing for the nonsinglet 



N + 1 



+ ^Cf + 0{N + 1) 



(2.15) 



(2.16) 



where Ca = 3, Cp = |, Tji = ^ for QCD with three colors. The corresponding splitting 
functions are simply found by noting that 



M[d{l~x)] = 1, 



M[l] 



N+l ■ 



(2.17) 



Positive powers of correspond to logarithmic derivatives of 6{l — x) which, when summed 
up, produce the nonlocal x propagation kernel of eq. (|2.3|) . If instead the expansion is 
truncated, the evolution equations are local. Specifically, using the splitting functions 
obtained according to eq. ( |2.17|) from the anomalous dimensions eq. ( p.l5|) in the evolution 
equations eq. (|2.3|) , and then differentiating with respect to ^, we get 



d 



d 



2 



-qCa 



2 
3 

^TRUf 



d_ 



+ 



2C^ 




2Cf 




(2.18) 



Cf_ 



' d_ 



+ 2 



9(^,0, 
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where we have for convenience defined 



(2.19) 



The AltareUi-Parisi equations thus simply reduce to two-dimensional partial differential 
equations with constant coefficients, thereby proving the local nature of the evolution. 

The singlet evolution equations eq. ( |2.18| ) can be solved by diagonalizing the matrix of 
anomalous dimensions eq. (p.l5| ): the eigenvectors will then satisfy (decoupled) equations 
of the form eq. (|2.7|) , with anomalous dimensions given by the corresponding eigenvalues. 
To order A^'^ the eigenvalues are 



A- 



N 



11 , 
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CA+'iTnUf-^TnUf +0(iV); X_ = -^TnUf + 0{N) 



corresponding to the eigenvectors v± = {Q±tG±) given by 

Q+ = l^nfNG+ + 0(iV2) Q_ = -gAG_ + 0{N); 



(2.20) 



(2.21) 



the singlet quark and gluon distributions are then given by Q = + Q- and G = 
+ G-. Notice that, according to standard perturbation theory, when the eigenvalues 

are determined up to next-to- leading order (in N) , only the leading nontrivial term ought 

to be kept in the expression for the eigenvector. 

Transforming to x space with the help of eq. ( |2.17| ) we thus get 



d 

d 
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+ 6. 



~d_ 



G+(e,0 = o, 
^-(^,0 = 0, 
fl 9(^,0 = 0, 



(2.22) 



with 



7 



12 

Jo'' 



3/3o 



11 + H!V 
±± -r 27 



16n/ 

Wo 



Wo 



(2.23) 



The eigenvector conditions eq. ( p.21|) become 

nf dG+{iX). 
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G_(e,C). 



^2.24) 



The eigenvector condition for can be conveniently rewritten by differentiating both 



sides with respect to C, and using the expression eq. (|2.22|) for -q^qq'- the term proportional 

dG 

to can then be neglected since it only gives a subleading correction to the eigenvector, 
and we get 
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(2.25) 
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Only the "large" eigenvalue A_(_ is singular at small A^, hence, at small x and large Q^, 
the "small" eigenvalue may be neglected. The singlet quark and gluon distributions are 
then no longer independent, but rather related by the large eigenvector condition eq. ( |2.24|) . 
If we had kept only singular terms in the anomalous dimensions, the small eigenvalue would 
have vanished, the large eigenvector would have coincided with the gluon, and the singlet 
quark distribution would have vanished. This is also apparent by direct inspection of the 
small X evolution equation ( p.22| ), which would have coincided with eq. ( p.l8| ) with all 



terms proportional to ^ on the right hand side neglected. When the constant terms are 
retained, the large eigenvector contains a mixture of gluons and singlet quarks, obeying 
the evolution equation eq. ( |2.22|) for G+ (and Q+). This has the same form as the gluon 
equation in eq. ( |2.18| ), but with the inhomogeneous term neglected, and a slightly different 
value of the coefficient of The singlet quark contribution Qj^ is however subleading, 
being determined in terms of the gluon contribution according to eq. (|2.25| ), which 
coincides with the singlet quark equation in eq. ( p.22|) . 

It thus appears that the behaviour of the singlet and nonsinglet components of F2 at 
small X and large are entirely determined by those of G'+(^, C) ^ind g(^, () respectively 
which, in turn, are determined by the evolution equations eq. ( p.22| ) . These are recognized 
as a two-dimensional wave equations, i.e. two-dimensional Klein-Gordon equations written 
in light-cone coordinates Q = x ± t with imaginary mass. This immediately implies 
several general properties of their solution: 

(i) The equations are essentially symmetrical in ^ and so G+{^X) ^-nd q{^X) evolve 
('propagate') equally in both ^ and ( (i.e. in x and Q^), up to the (small) asymmetry 
induced by the 'damping' term proportional to 6. Any further asymmetry in ^ and ( 
must thus come from the boundary conditions. 

(ii) The propagation is 'timelike', into the forward 'light-cone' at the origin (^, Q = (0, 0), 
along the 'characteristics' ^ = constant and ( = constant (see fig. 2). 

(iii) At a given point (^, (), G+(^, C) ^ind q{^, () depend only on their respective boundary 
conditions contained within the backward light cone formed by the two characteristics 
through (^, C)- 

(iv) Because the equations are linear, contributions to G+{^X) or q{^X) from different 
parts of the boundaries are simply added together. 

(v) Since the 'mass' terms are negative, the propagation is 'tachyonic'; this means that 
both G+{$,X) and q{CX) are unstable, growing exponentially rather than oscillating. 

(vi) Since 6 and d are both positive, the damping terms ensure that, at fixed ^, G'_|_(^, () 
and g(^, C) eventually fall with increasing (. 

It is also straightforward to obtain the general solution to the evolution equations: 
they are simple examples of the characteristic Goursat problem, in which the solution is 
entirely determined by the knowledge of boundary conditions along two characteristics, 
and can be written explicitly in terms of a particular solution to the equation. The latter 
is easily found by observing that, setting z = 27v^, the first and third of eqns. ( |2.22| ) 
coincide with the Bessel equation, the appropriate solution of which is the Bessel function 

^o(^)-f;%n^ ~ 7^^^1 + 0(1))- (2-26) 
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The general solutions are thus found to be 



G+it = /o(27v/^)e-^+^G+(0, 0) + Io{2^Vif^)e-'+^ -^G+i^ , 0) 

G'_(^,C) = e-^-^G_(^,0), 

= /o(27 VlC)e-^'^Q(0, 0) + d^'/o (27V(e-e')C)e-^^^?(e', 0) 

+ 1^ dC'lo (27V^(C-C'))e^(^'-^) (^?(0, CO + 5~q(0, C')) • 

(2.27) 

These solutions display explicitly the symmetric nature of the evolution. It is interesting 
to compare them with the standard solution to the Altarelli-Parisi equations at larger x, 
where boundary conditions are imposed at a scale to for x' such that x < x' < 1. Since the 
structure function vanishes kinematically at a; = 1 the boundary condition on the lower 
boundary is trivial, and evolution takes place from the initial parton distributions assigned 
at a given to forwards in t. In the present case, instead, the two boundaries are treated 
symmetrically, and evolution takes place as much with respect to x as it does with respect 
to t. 

The asymptotic behaviour of G_|_(^,C) or q{CX) small ^ and large C, i-e. far away 
from the boundary, will in general depend on the form of the boundary conditions. Due 
to the linearity of the equation, we may consider each boundary separately. Contributions 
from each boundary are generated by fluctuations of the functions or q on that bound- 
ary. If these fluctuations are sufficiently well localized close to the origin, then far from 
the boundary we can use a multipole expansion, expanding the argument of the Bessel 

functions in powers of ^ (i.e., the distance from the boundary over the spread of the 

a/ 

source) for the left boundary and ^ for the lower boundary. All the contributions from 
higher moments of the boundary fluctuations are then seen to be suppressed by powers of 
the light-cone distance from the origin 



a = VeC, (2.28) 

while the leading contribution is simply given by the strength of the source at the origin, 
and its asymptotic behaviour is determined by that of the Bessel function (|2.26| ): 

G4p, o) ^ Ar-gL= exp {2ja - (5) } (l + 0(i)), 

where we have introduced the hyperbolic coordinate orthogonal to a, namely 

P = VUC (2.30) 



(2.29) 
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(see fig. 2). 

The origin of double scaling is now clear: because of the isotropy of the evolution 
equation, its solution asymptotically only depends on the scaling variable cr, i.e., its level 
curves are hyperbolae in the (^,C) plane (see fig. 2), and do not depend on the direction 
in which the propagation occurs [p scaling). Furthermore, because the solutions ( |2.29| ) 
are asymptotically independent of the boundary conditions, the dependence on a is given 
by a universal rise, stronger than any power of ^ but weaker than any inverse power of 
X. The universal form of this rise reflects the underlying dynamical mechanism which 
generates it, namely, in the singlet case the collinear singularity in the triple gluon vertex 
which is responsible for the singularity in the anomalous dimension (eq. (|2.15 )) whose 



strength determines the coefficient 7^ (eq. ( |2.23| )), and in the nonsinglet the corresponding 
singularity in 7^^ due to collinear gluon bremsstrahlung, whose strength determines the 
coefficient 7^. In an abelian theory the former singularity clearly vanishes, so the singlet 
distributions would no longer grow so strongly, behaving instead more like the nonsinglet. 

Using eq. ( p.25| ) to determine Q, and the expression eq. ( |2.5| ) of F2, and neglecting 
the nonsinglet contribution, since this is suppressed by a factor of approximately e~^, it is 
easy to show that the asymptotic behaviour of -F|' (or F2) will also have a double scaling 
form, namely: 

- ^-G+(a,p)(l + 0(i)+0(i)), (2.31) 

with given by eq.( |2.29| ). This behaviour holds up to corrections of order ^, from 
the subasymptotic form of the Bessel function and the boundary corrections, up to terms 
of relative order ^ in the exponent, from higher order contributions to the small 

expansion eq. ( |2.15| ) of the anomalous dimensions, and up to corrections of order i, from 
higher order contributions to the eigenvector equation. It thus holds in the limit a 00 
along any curve such that also p ^ 00, such as for example the curve ^ oc C^^"^ with e > 0: 
that is, far from the boundaries, and provided the increase of In - is more rapid than that 
of Int. 

All of this, however, hinges on the assumption that the solution eq. ( |2.27| ) to the 
wave equation may be treated in the multipole expansion, i.e., that the fluctuations on 
the boundaries fall off away ^from the origin. If this does not happen, then the multipole 
expansion is not valid, the asymptotic behaviour eq. ( |2.2y| ) does not hold, and the chain 
of arguments leading to the scaling form of F2 eq. ( |2. 3 1| ) breaks down. Specifically, assume 
for example that the boundary conditions at ^ = are exponentially rising functions of 
^: G'_|_(^, 0) ~ expA^ = x~^, ?(^7 0) ~ expA^ = x"'^. The boundary integral may then 
be evaluated by the saddle point method, is dominated by a nontrivial saddle-point, and 
gives the asymptotic behaviour 

G+{a,p) ~ Ar'exp{Aap+(:i^-5+)(^)}(l + 0(i) + 0(i)), 

(2 32) 

q{a,p) ~ Ar'exp{Aap+(£-5)(^)}(l + 0(i)+0(i)). 

(T — >00 /\ t-' 

The strong growth on the boundary is thus preserved by the evolution. Since the evolution 
equations are linear, this behaviour should be added to the dynamically generated con- 
tributions (|2.29| ). Since it is powerlike, the rise of the boundary condition will eventually 
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dominate the dynamically generated rise eq. ( 2.29 ) when ^ is large enough. However the 
universal behaviour eq. (|2.29|) is still dominant when ( is large enough that the nontrivial 
saddle point leading to eq. (|2.32| ) is no longer dominant, i.e. whenever p ^ p ^ ^ r&- 
spectively. Similar contributions to ( pj.32D would arise from exponentially rising functions 
of C on the lower boundaries. 

In practice, we may choose xq close to the turning point of the evolution, so that 
the lower boundary condition is reasonably flat. Furthermore Regge theory suggests that 
the left hand boundary condition is given by F2 ~ x^~°'^^\ where a(0) is the intercept of 
the appropriate Regge trajectory. In the singlet channel, this is the pomeron trajectory, 
with ap(0) ~ 1.08, so the left hand boundary condition for G+ should also be soft, and 
the asymptotic behaviour of F2 and F2 given by ( p.31| ), at least for p ^ 7/(a;p(0) — 1) ^ 
15. However in the nonsinglet channel, the appropriate trajectory is that of the p, and 
p(0) ~ ^. So the left hand boundary condition for the nonsinglet, q ~ x""'', is hard, and 



a 



asymptotically ( |2.2y| ) will be dominated by ( |2.32| ), whence, for p ^ 1 



^2 



M" exp {(A - l)ap + (f - )} (1 + 0(i) + 0(i)), 



(2.33) 



with A ^ ^. Recent data from NMC [|T2| seem to be in good qualitative agreement with 
this prediction, but deuteron data in the kinematic region explored at HERA would be 
necessary to confirm it more precisely. 

It was suggested some time ago |T^ that the singlet boundary condition at Qq might 
not be given by the intercept of the pomeron trajectory, but could rather rise very steeply 



as x~^^, with \l = 41n2^a, 



1 

2- 



This steep initial rise (sometimes called the 'hard 



pomeron') was supposed to incorporate. 



higher order perturbative effects at small x described by the BFKL equation [P 



in an admittedly rather heuristic manner, the 

With 



such a boundary condition, the double scaling rise (|2.31|) in Ff' would be masked by the 
stronger rise of as given by ( p.32|) , so 



^2 



Ar"exp{A^ap+(^-5)(^)}(l + 0(i) + 0(i)). 



(2.34) 



whenever p~7/Al — 2. This prediction is nonuniversal, in the sense that the precise slope 
of the rise cannot be predicted since A^ depends on as, and it is not known at which scale 
Us should be determined. Furthermore it is qualitatively different in form from the double 
scaling rise (2.31): in particular the rise at large ^, fixed C, is now no longer accompanied 



by a corresponding rise at large C, and fixed (though large) ^. We will come back on the 
issue of applicability of leading order computations and the relative importance of higher 
order corrections in the next section. First however we will see whether the HERA data 
support the universal double scaling prediction (|2.31| ) , or prefer the more phenomenological 
suggestion ( p.34| ). 

The scaling plots in fig. 1 display the measured values of F2, with the subasymptotic 
corrections in eq. ( |2.31| ) rescaled out, i.e. R'pF2 with 



R'f{(^^ p) = exp {5{a/p) + ^ ln7a + ln(p/7)). 



(2.35) 



10 



The scaling variables are computed with xq = 0.1, Qo = 1 GeV, A = 263 MeV, and S = S+ 
eq. ( |2.23| ) with Uf = 4. The data cover a wide span in p: for instance the data in fig. Ic 
have 1 Nevertheless they all fall on the same line, and display a slope which agrees 

very well with the predicted asymptotic value 27 = 2.4. If the boundary condition were 
hard, eq. ( |2.32|) shows that the leading behaviour would also be a linear rise of lnF2 in a, 
but now with a slope which is not universal (as it depends on A), and strongly p dependent. 
Hence, the data should not fall on a single line, and the agreement of the observed slope 
with the calculated value of 27 could only be a coincidence. The fact that the data display 
double scaling thus allows us to exclude the possibility of power-like boundary conditions 
to leading-order (or, as we will se in a moment, next-to-leading order) evolution at a very 
high confidence level 

The scaling plots in fig. 1 also show that the slope of the rise of In F2 is significantly 
smaller than the asymptotic one when a is not too large. This suggests that scaling 
violations may already be important here. Even more dramatic scaling violations are seen 
if one considers data at low Q^. Both effects are illustrated in fig. 3 (i), which displays F2 
after complete rescaling of the leading asymptotic behaviour, i.e. by a factor 



Rpicr, p) = exp ( - 270- + d{a / p) + \ In 70- + ln(p/7)). 



(2.36) 



It is apparent that the data display a systematic drop in a (the asymptotic double scaling 
line is approached from above), and the recent (albeit preliminary) data with low 
(2 ^ p ^ 5 and 1.3 ^ cr ^ 1.8) do not seem to scale at all. 

This leads us to consider scaling violations, the simplest of which appear when two- 
loop corrections are included in the Altarelli-Parisi equations. The leading singularities of 



the two loop anomalous dimensions |jT^ in MS scheme are 



1 

iV 



(|Cf 



-gCATRUf 

+ 0((Ar + l)-2). 



fCpTRUf ) 



+ 0(1), 



(2.37) 



Note that all the entries in the two loop singlet anomalous dimension are singular, and 
the nonsinglet is more singular than at one loop. Taking the inverse Mellin transform of 
( p.37|) , using for the nonsinglet the result In^ \~\ = 1/-^^, and including them in the 

Altarelli-Parisi equations ( |2.3| ) gives wave equations similar to (|2.18[ ), but with additional 
terms of 0{as) on the right hand side. Linearizing these two loop corrections, the wave 
equations (|2.22[) become (suppressing subasymptotic contributions) 



d 



+ 5 



G'_(^,C) = e_a,(to)e-^G+(^,C), 



^2.38) 



Q(^,c) = e«,(to)e-^ / dei^-aqi^'x) 
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with now, in place of (|2.11|) , 



C^ln 



(2.39) 



with as evaluated to two loops, and 

103 ^(3i 



e+=(^n, + 3-^)/7r/3o, 6_ = 



26n/ 
37r/3o' 
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(2.40) 



(3i = 102 — ^rif being the two loop coefficient of the /9-function. It is not difficult to 
show that the three new parameters ( 2.40 ) are in fact all independent of the choice of 
factorization scheme. The eigenvector conditions ( p.21| ),( p.24|) and (|2.25| ) are unchanged. 

The general solutions of the three equations ( |2.38| ) are the same as those of the leading 
order equations, ( ^.27| ), but each with an additional contribution on the right hand side: 



Po Jo Jo 

e.^ /'<e^-(^'-^)-^'G+(r,C'), 
Po Jo 



(2.41) 



An 
Po Jo 



JO 



' dC'dedC'Io (27 Vii-OiC-O) e 



S+{c'-C)-C'(c' 



respectively. Two loop corrections to the leading asymptotic behaviour are then found by 
substituting the leading behaviour into ( |2.41j ), and evaluating the asymptotic form of the 
integrals. Clearly the form of the correction will then depend on the form of the leading 
behaviour. For soft boundary conditions the double scaling behaviours ( |2.29| ) are corrected 
by extra factors 



G+(p,a) 



Af- 



1 



.^oc - • y4F^ {2^^ - (f ) ) [1 - ^+ i^sito) - «.(t)) (1 + 0(i)) . 
q{p, a) ~ Ar^=L= exp hja -~5{^)\\l- e(a,(to) - a.(t)) 41 (l + 0(i)), 



while the double scaling behaviour of F2, ( |2.31| ), is corrected by a similar factor: 



(2.42) 



(T^oo 162 p y/ATTja 



1 - {e+ {asito) - as{t)) - ^a.(t)) £ (l + 0(i) + 0(i)), 



^2.43) 



the extra term in the square brackets coming from the subleading contribution to G- . For 
hard boundary conditions, the asymptotic behaviours ( p.32|) , ( |2.33| ), are corrected by the 
same factors ( p.42| ), ( |2.43| ) respectively, but with p/7, p/7 replaced by A, A respectively: 
the corrections are thus then x independent. 
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Two loop corrections will be most important for soft boundary conditions, and thus 
in the singlet channel, as the leading correction to the double scaling seen at HERA. To 
see their effect, the data are replotted [in fig. 3 (ii)] with a new rescaling function. 



d(2), 
tip 



(a, p) = Rpia, p) [l - a,(to) (e+ - (e+ + ^)e-^'') f ] (2.44) 



where Rp is the leading order rescaling ( |2.36|) . It can be seen from the plots that the 
effect of the two loop correction is moderate (except at very low Q^) but significant in the 
range of the present data: it increases the starting scale from Qq ^ 1 GeV at leading order 
to around Qq ~ 1.5 GeV, reduces the slope of the a-plot by about 10%, and decreases 
the rise in the p-plot at large p and low Q-^. Thus most of the as yet observed scaling 
violations can be accounted for by the two loop correction, and conversely the effect of this 
correction can be clearly seen in the data. These results have been confirmed by numerical 



calculation using the full one and two loop anomalous dimensions [ 1^-19 ]■ 

The residual rise at large p in the low data fl^, if it turns out to be statistically 
significant, could be due to many different (and possibly competing) nonperturbative ef- 
fects: a small rise in the (nonperturbative) boundary condition (even x~^'^^ has observable 
effects at Qq = 1.5GeV), nonperturbative effects due to the opening of the charm thresh- 
old, conventional higher twist effects, or even more novel higher twist effects such as parton 
recombination. However it could also be due to higher loop singularities. It is particularly 
important to consider these, since it is necessary, in the light of the above discussion about 
the 'hard pomeron' boundary condition above, to understand why they do not in fact spoil 
double scaling. Such an undertaking is possible since the precise form of the leading (and 
some of the sub-leading) singularities are known: we will now explain in some detail how 
their effects may be properly included. 



3. Perturbation Theory at Small x 

In the previous section we have seen how double scaling appears as a generic feature 
of the solution to the LO or NLO Altarelli-Parisi equations, in the limit as p and a grow 
large, i.e. as both In ^ and Int grow, provided the former grows faster than the latter. This 
prediction thus defines a "double scaling limit" of QCD, intermediate between the Regge 
limit (small x at fixed t), and the Bjorken limit (large t at fixed x). In the Regge limit 
perturbation theory fails and we are unable to calculate either the x or the t dependence 
of parton distributions. In the Bjorken limit perturbation theory holds and predicts the 
scale dependence of parton distributions, viewed as functions of x, but the x dependence 
itself depends on an uncalculable initial condition. In the double scaling limit, the x and 
t dependence is universal and only depends on a single overall normalization. 

We can thus divide the (^,C) plane in various regions (see fig. 2): for low ( (say, 
C < with a suitable choice the origin of coordinates) perturbation theory breaks down, 
whereas the Altarelli-Parisi equations hold for positive values of C, and become more 
and more accurate as ( increases. On the other hand if ^ is also sufficiently large, the 
anomalous dimensions may be expanded around their leading singularities eq. ( |2.15D so the 



Altarelli-Parisi equations take the small x form eq. (|2.18| ). However if ^ keeps increasing. 
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perturbation theory eventually breaks down because higher twist corrections, necessary to 
ensure unitarity, must become get more important. Physically, one can understand this by 
noting that x can be interpreted as the momentum fraction carried by individual partons. 
When ^ is very large momentum is shared between an increasingly large number of partons; 
this corresponds to an increase of the parton density which cannot continue indefinitely 
and should eventually stop when partons start to recombine with each other This 
is expected to happen in the region x ~ Xr exp{—as{to)^ /asit)"^), where Xr cannot be 
reliably computed but could be of order 10~^ at scales of a few GeV^ If instead ( 

keeps increasing we eventually get to a region where effects from the lower boundary of 
perturbative small x evolution are important: these reflect the shape of parton distributions 
propagated down from large values of x, and thus we get back to Bjorken scaling, with 
uncalculable dependence on x. 

Even in the remaining double scaling region, however, double scaling actually holds 
only with sufficiently soft boundary conditions; besides, we have derived it only in a LO or 
NLO calculation. The two issues are actually closely related, because of general arguments 
suggesting that when higher order corrections are included a hard power-like rise of the 
singlet parton distributions may result. The underlying logic is the following: the leading 
logarithmic gluon-induced contributions to the deep-inelastic scattering cross section to all 
orders in ag log ^ may be summed up by the solution of an equation satisfied by the gluon 
distribution (BFKL equation) |T^. This equation is not consistent with the renormaliza- 
tion group, in that it does not sum logs of and thus it does not include evolution in Q^: 
in fact, it is derived at fixed coupling a^. It therefore also does not separate leading twist 



(3.1) 



from higher twist contributions. However, it is possible to extract the leading twist 
contribution to the solution of this equation in Mellin space: as oo the solution has 

the form 

the extra term in the square brackets being higher twist. The leading twist piece is now 
the same as a solution to a (fixed coupling) renormalization group equation eq. (p77|), with 
7(A^, as) identified as the anomalous dimension. This anomalous dimension is determined 
as the inverse of the function x(^) = 2'i/'(l) — i(j{x) — — x) (where 'iIj{x) is the Euler 
function) : 

x[7(iV,as)] = — , (3.2) 



a. 



where a. 



The solution to (|3.2|) , the Lipatov anomalous dimension 71,, thus turns out to be a 



function of -^: 



,,(../iv)=i:.i"(|)' 



(3.3) 



k = l 



with coefficients jj^^ determined uniquely by eq. ( p.2| ). All the coefficients turn out to be 

positive, save 7}^^'' and 7^^'' which vanish. Now, inverse powers of are the Mellin 
transforms of logs of - : 



M 



"1 1 

x{k-l)\ 



In 



fc-i 



X 



1 



(3.4) 
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One may then argue pO] that, since the leading logs of - are summed by the BFKL 



X 

equation which leads to eq. (|3.1j), the coefficients 7|^'' in eq. (|3.3| ) must give the coefficient 
of the most singular term in expansion in powers of N of the fc-th order contribution 
to the ordinary anomalous dimension (leading singularity): these are by definition the 
leading logs in ^ since in x space they correspond to the contributions with the largest 
number of logs at each perturbative order. This result may actually be proven rigorously 



by means of suitable factorization theorems [|211] : the expansion of the Lipatov anomalous 



dimension eq. ( p.3| ) gives the coefficients of the leading singularities to all orders in ctg in 
the gluon anomalous dimension (due to the fact that the BFKL equation describes 
gluon propagation and emission). Notice that these coefficients are factorization scheme 
independent. 

This result has several important consequences for our discussion. First, it shows 
explicitly that the double logarithms, of the form -as{cts In'^ -)"^~^, which one might expect 



naively to arise in a perturbative expansion of splitting functions |2^], are reduced in 
the gluon sector to single logarithms of the form ^as(as In i)"^~^, because many of the 
singularities cancel systematically. This cancellation actually occurs in the whole singlet 
sector (though not in the nonsinglet [^). This is in accordance with the calculation 
of the two loop singularities ( |2.37| ) (which in fact exhibit a further accidental cancellation 
in the singlet sector, which also occurs at three and five loops). 

Despite this remarkable cancellation, at higher orders in ctg the singularity in the 
anomalous dimension is still growing strongly, albeit not quite so fast as one might have 
naively expected. Thus at small enough x the enhancement due to the extra powers of In ^ 
in the corresponding splitting function may offset the suppression due to the extra powers 
of Og, so that the inclusion of higher order corrections may be required in order to obtain 
accurate results. Furthermore, the function x(x) has a symmetric minimum at x = 1/2, 
implying that the anomalous dimension 'jL^cts/N) has a square-root branch point there 
(there are also other branch points, and in fact the structure of ^l{ois/N) in the complex 



plane is quite complicated |]2^). The value of its argument such that 'fLicts/N) = 1/2 is 
4 In 2, so the branch point is at 

AL(a)^41n2f^Y (3.5) 



TT 

This implies that Mellin-space parton distributions also have a singularity at = Al , and 
correspondingly, in x space they should increase as x — > as x^^. This corresponds to a 
rather strong rise for realistic values of ctg, and would certainly spoil the observed double 
scaling behaviour if it occurred already at presently attainable values of x. However, it 
is not clear that this simple argument is consistent with the renormalization group: for 
example, the result of using a power like behaviour of this form as a boundary condition 
for LO or NLO perturbative evolution depends on the scale Qo at which the boundary 
is set. Moreover, it is not clear how this behaviour should be matched to conventional 
perturbative evolution at larger x, nor indeed in which region of the x-t plane it sould 
become important. 

We need thus to keep into account higher order singularities in in a way consistent 
with the renormalization group. This can be done |]2^ by reorganizing the perturbative 



expansion of the anomalous dimensions, or, equivalently, the splitting functions used in 
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the Altarelli-Parisi equations. To understand how this works, it is convenient to classify 
the contributions to anomalous dimensions by expanding the anomalous dimension used in 
renormalization group equations eq. (|2.7|) in powers of ttg, and then each order in powers 
of N (see fig. 4): 

oo m oo ^ m \ 

g7(iV, «) = E E iV-^ = E ( E ^nN-- + l^N^ ) , (3.6) 

m=l n= — oo m=l ^ n=l ^ 

where the numerical coefficients A"^ are given by eq. (|3.3|) , and in the last step we have 



separated out the regular part of the anomalous dimension Using eqs. (|2.17|) , (|3.4|) 

this is seen to corresponds to expanding the associated splitting functions as 

oo / 1 1 '""^ i \ 

= 5: („, -( i 5; + (..)), (3.7) 

m=l ^ 71=1 

where P^™-) (x) are regular as x — > 0. 

Solving renormalization group equations sums all leading logs of the scale which ap- 



pears in the equation: for instance, upon solving eq. (2/7), the anomalous dimension gets 



exponentiated according to eq. (|2.1UD. At LO only the term with m = 1 is included in the 



anomalous dimension eq. ( |3.6| ), which is thus linear in ctg, so this amounts to summing 
up all contributions to the deep-inelastic cross sections where each extra power of ag is 
accompanied by a power of InQ^. In fact, because the LO anomalous dimension has a - 
singularity (which leads to a factor of In - in the cross section upon integration) some logs 
of - are also summed, but In - is not considered leading, in that factors of as may or may 
not be accompanied by In ^. Thus, if the LO in the expansion eq. ( p.6|) of the anomalous 
dimension is used, all logs of the form 

aP(lng2)«(lni^ (3.8) 

with q = p, and < r < p are summed up. At NLO the anomalous dimension includes 
both linear and quadratic terms in ctg, and thus all logs eq. ( |3.8| ) with q < p < 2q are 
summed (there are at most twice as many powers of as as powers of InQ^). In standard 
perturbative computations the NLO solution is however then linearized, by expanding the 
exponential of the NLO term of the anomalous dimensions in powers of as and retaining 
only the leading nontrivial term (for instance, the NLO asymptotic correction eq. ( |2.43| ) 
was derived in this way). This means that only terms with q < p < q + 1, i.e. p = q + 1 
are included at NLO. Furthermore, there is an extra power of In ^ in the NLO anomalous 
dimension, so < r < p as at LO. At NNLO yet an extra power of as per power of InQ^ 
is allowed, and so forth. We will refer to this as the "large x" expansion. 

This is not the only way to organize the perturbative expansion, however. We might 
instead want for instance to consider In ^ as leading: this would be appropriate at very 
small X. Then, all terms where each extra power of as is accompanied by a power of 
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In ^ should be included at LO. These are the leading singularities to all orders in ag (see 
fig. 4b) : it is thus convenient to reorganize the expansion of the anomalous dimension as 



oo oo / \ 

l^l{N.a)=J2a-\ Y: ^r-^-M^j ), (3.9) 



m=l n=2 — m 

which corresponds to the splitting function 

oo / ^ oo / n—1 1 



m=l n=l 
m-2 



+ 5:^-— (a.(t))-|^5(l-x)]. 



(3.10) 



5=0 

Subsequent orders are still labelled by the index m of the outer sum. The LO is the sum 



of leading singularities, and sums all logs eq. (3^) with r — p and 1 < Q < p. (terms with 
q — Q are not included because at least one power of InQ^ is produced by integration of 
the renormalization group equation) . In NLO terms with an extra overall factor of as are 
included in the anomalous dimension, so the solution contains terms with r < p < r + q. 
Upon linearization, only terms with p = r + 1 are kept, while still 1 < q < p. There is then 
complete symmetry between this "small a;" expansion and the large x expansion, with the 
roles of the two logs interchanged. 

The small x approximation to LO evolution discussed in the previous section cor- 
responds to only retaining the most singular terms in the LO anomalous dimension: it 
thus corresponds to taking the intersection of the small x and large x leading order terms, 
i.e. the pivotal term with m = n = 1 in eq. ( ^.7] ) or ( |3.1UD , which sums all logs with 
p = q = r. The symmetry of double scaling refiects this double logarithmic approxima- 
tion. The terms contributing to the coefficient 5 eq. ( |2.23D are large x corrections (i.e 
corresponding to r < p = q) and so forth. The fact that double scaling is observed indi- 
cates that the HERA data are taken in a region where the two logs start being equally 
important. 

This suggests that in this region the most convenient way of organizing the perturba- 
tive expansion is one where the two logs are treated symmetrically at each order ("double 
leading" expansion). To do this, the anomalous dimensions are expanded as (see fig. 4c) 



a 

£7(iV,«) 

ZTT 



oo ^ oo / \ \ 



In LO (i.e. when m = 1) all terms with l<q'<p, 0<r<p, and 1 < p < q ^ r are 
summed. If all cross terms (i.e. those containing a product of a contribution to 7^'* times 
a singular contribution) are linearized then the solution includes all contributions where 
each extra power of as is accompanied by a log of either ^ or or both. In NLO an 
overall extra power of as will be allowed, and so forth. 
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Of course, a variety of other expansions which interpolate between these could be 
constructed. The crucial point here however is that all these expansions are consistent 
with the renormalization group: in each case we may define 



^7(A^,«) = «"'7m(iV,«); (3.12) 

each term in this expansion is then of order as compared to the previous one and, in par- 
ticular, a change of scale at /c-th order may be compensated by adjusting the k + 1-th order 
terms. The expansion eq. (|3.12D may thus be treated using the standard machinery used 
to perform NLO or higher order QCD computations. However, if we choose an expansion 
which is appropriate at small x (say, the extreme small x one ( p.9|) ), then 7^+1 (A^, a) is 
of the same order as 'jmiN, a) since they both sum up the relevant logarithms, hence the 
m-fl-th order contribution to 7 is genuinely of order ctg as compared to the m-th order one. 
This is not achieved by an ad hoc "resummation" of a particular class of contributions, 
but simply by organizing the perturbative renormalization group in a different, equally 
consistent way. 

There is a very important issue which must still be addressed however: namely, in 
the large x expansion eq. (|3.6|) the series in n which defines the anomalous dimension 
at each perturbative order (in m) of course converges — in fact, we defined this series 
by expanding out an expression given as a function of A^. The summation of leading 
singularities eq. ( |3.3| ), however, does not converge for all A^: as mentioned earlier, 7l(A^, cts) 
has a branch-point, implying that the series has a finite radius of convergence, and diverges 
when < 41nQ; (the other two branch points on the first sheet being also inside this 



circle). In fact it is not even (Borel) resummable for ReA^ < 4 In a, since the integral 
over the Borel transformed series (which has infinite radius of convergence) diverges at the 
upper limit. This seems to pose an insurmountable problem for the perturbative approach 
to small X evolution: the series which defines the leading coefficient in expansion (|3.9| ) of 
the anomalous dimension, which was supposed to be useful for small A^, is actually only 
well defined when A^ is sufficiently large. This apparent inconsistency seems to have led 
many to the conclusion that conventional perturbative evolution breaks down at small x. 

The dilemma is resolved by the observation that the physically relevant quantity, 
namely the splitting function Pgg{x), is instead given at leading order in the expansion 
( p.lOp by the series 

^^E(-^^)°° '"-nr' - '-^HKmni), (3.13) 



n=l 



X 



where X{as) is given by eq. (^ 



B(-) = E7^-"-' (3-") 



n=l 



a^ = {A^y^{4ln2CAM-^. (3.15) 
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Because the series Yl'^=i O'nV'^ has radius of convergence one, the radius of convergence of 
B{u) is infinite, and thus the series which defines the sphtting function eq. (|3.13| ) converges 
uniformly on any finite intervals of x and t which exclude a; = 0: the only reason for the 
bad behaviour of the series ( p.3|) is that when transforming to Mellin space one attempts 
to integrate all the way down to x = 0, and this is not possible for singlet distributions 
because the total number of partons diverges there. If instead the parton distributions are 
evolved using the Altarelli-Parisi equations ( p.3| ), the splitting functions are only required 
over the physically accessible region x > Xmin, and no convergence problems arise. Indeed, 
because the series for the splitting function is convergent it is only necessary, when working 
in a physically accessible region to a certain level of accuracy, to retain a finite number of 
terms. Since then the only singularities in the Lipatov anomalous dimension are poles at 
the origin, it follows that the cuts which arise to all orders are strictly unphysical. 

Similar considerations will apply to the subsequent orders in the expansion of the 
anomalous dimension eq. ( p.9|) : if the coefficients in the expansion of the anomalous di- 
mension have a nonvanishing radius of convergence, then the corresponding splitting func- 
tions will converge for all a: > 0. This will be true provided there is no singularity (or 
accumulation of singularities) in the anomalous dimensions at = cxd or, equivalently, at 
as = 0. This is a standard assumption in perturbative QCD (at least order by order in 
perturbation theory), so it is natural to conjecture that this will be true for all values of 
m in eq. (W^)- 

So far we considered the behaviour of the anomalous dimension. We already know 
that at LO only and 7^^ display a singularity at N = 0. In fact, it may be proven [^| 
that the leading singularities in •jgg and 7^5 are related to all orders, according to the 
so-called color-charge relation 



7S9 



Cf 
Ca 



I99 + O 



a. 



iV 



(3.16) 



while the coefficients of the leading singularities in the quark anomalous dimension 7^^ 
and 7qg vanish to all orders. The coefficients of the NL singularities (i.e. the NLO in the 
small X expansion eq. ( |3.9|) ) have been calculated recently [^: beyond the (nonsingular) 



lowest order in [given in eq. (|2.15| )] they also satisfy a color charge relation, namely 



Iqq 



Cf 
Ca 



a,. 



OLs 



(3.17) 



The series expansions for these NLO anomalous dimensions have the same radius of con- 
vergence as that of the Lipatov anomalous dimension, though the form of the branch point 
singularity is now scheme dependent. The NLO coefficients in the expansion of the gluon 
anomalous dimensions are still unknown. 

All the NLO coefficients, being subleading, are scheme dependent: given their expres- 
sion in the parton scheme, where coefficient functions satisfy eq. ( |2.2D , their expressions 
in any other factorization scheme such as the MS scheme may be determined, along with 
the corresponding coefficient functions. In fact, the freedom of choosing a factorization 
scheme turns out p8|| to be wider in all expansions other than the standard large x ex- 
pansion eq. (|3.6| ): this is due to the fact that the normalization of the gluon beyond order 
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as may be modified by change of sdieme. It follows that, besides the usual freedom to 
perform a scheme change which modifies the F2 coefficient functions, there is now also the 
possibility of performing a scheme change which does not affect the F2 coefficient func- 
tions, but changes the definition of the gluon distribution (while leaving all LO anomalous 
dimensions unaffected, as a scheme change ought to). Such a redefinition, in particular, 
affects the quark anomalous dimensions, which only start at NLO in the small x expansion: 



for instance, it can be used [29| to remove a singularity which this anomalous dimension 



has at NLO in the parton scheme [|2^ at the Lipatov point N — Ai,(as)- The nonsingular 



anomalous dimensions thus obtained has been argued |29] to have a more direct physi- 
cal interpretation in that they correspond to specifying the initial parton distributions in 
accordance with Regge theory. It is also possible to find a scheme in which the quark 
anomalous dimensions vanish [pO| . 

This freedom, however, could in principle be completely pinned down by requiring 
that the evolution equations conserve momentum, which implies the constraints 

7gs(l,«s) +733(1, "s) = 0, 7^g(l,as) +73g(l'«s) = 0, (3.18) 

in analogy to what happens at large x, where the gluon normalization (which in that case 
is just a constant) is fixed by the same requirement. These constraints, which are necessary 
consistency conditions that must be satisfied order by order in the expansion of anomalous 
dimensions, cannot be satisfied at LO in the small x expansion, where however they do not 
apply because the the gluon decouples from F2 (which therefore does not evolve) . At NLO 
they fix uniquely the NLO coefficients in the small x expansion eq. ( p.9|) of in terms 



of the LO coefficients of •ygg and the NLO coefficients of [T^] : substituting eq. (|3.9|) in 
eq. ( p.l8|) implies immediately 



(K-^'U + (^n+')33 + iKXl)^9 = 0. (3.19) 



It can be proven that the freedom of choosing the gluon normalization by a change of 



scheme is in one-to-one correspondence with the values of the NLO coefficients in "fgg: of 
course, whatever the choice of gluon normalization there exists a set of values of {A'^'^^)gg 
which satisfies eq. ( |3.18| ); more interestingly, whatever the coefficients {A'^'^^)gg turn out 
to be when they will be calculated, there exist a choice of normalization that will enforce 
eq. ( |3.18| )." At present we have to content ourselves with parametrizing our ignorance of 
these coefficients by the choice of gluon normalization, and estimating the corresponding 
uncertainty by varying this normalization. 

In the nonsinglet channel the inclusion of the leading singularities is more problematic. 
The singularity at = — 1 in the anomalous dimension 7^^ is known to be stronger 



than that at A = of the singlet anomalous dimensions, i.e. to be double-logarithmic: 

2m-l 



"^7.T(^.«)= E"'" E ^niN + l)-''. (3.20) 



^ Strictly speaking this is only possible when n > 2 in eq. ( |3.18| ): the scheme-independent 
0{as) terms violate momentum conservation in the small x expansion. This violation is of course 
asymptotically subleading as x ^ 0, when this expansion is supposed to hold. Exact momentum 
conservation can be obtained order by order in the double leading expansion eq. 
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This is in agreement with the two loop result ( |2.37| ), and is indeed the generic expected 
behaviour of anomalous dimensions in the small limit. The behaviour of the nonsinglet 
distributions which is obtained by including these singularities to all orders in the coupling 
h as been computed |]2^, and leads to a power-like growth: qns ~ x~'^^^ with Akl 



y^Sa/Sn, with an almost identical result for charge-conjugation odd distributions. This 
is in fact very similar to the growth predicted by Regge theory and discussed in the 
previous section. However, it is obtained with fixed coupling and is not yet consistent 
with the renormalization group. Since the corresponding factorization theorems are as yet 
unproven, it is still unclear whether it will be possible to treat this case in the same way as 
the singlet one, making it consistent with the renormalization group by using the results 
of ref. to derive the leading term (namely the coefficients ^2m-i IHI]) ^^e small x 



expansion of the nonsinglet splitting function which can then be used in the Altarelli-Parisi 
equation in the same way as in the singlet sector. What is clear is that the effect of such 
a procedure would be relatively undramatic, given that the boundary condition is already 
hard, with A ~ Akl- It follows that the singlet contribution to F2 still dominates the 
nonsinglet at small x, and thus for the rest of the discussion we ignore it and consider only 
the singlet. 

We can now finally study the small x behaviour of singlet parton distributions: using 
evolution equations which include logarithmic effects in In ^ to all orders we will be able 
to assess in which region the simple LO double scaling predictions are reproduced. Even 
though in the HERA region the double leading expansion is probably more appropriate, let 
us consider the extreme small x expansion eq. ( |3.9| ), both because it is somewhat simpler, 
and because it allows us to test the stability of double scaling in an extreme case: if the 
LO correction in this expansion does not spoil scaling, no other corrections will since all 
the most singular corrections are included, and they all have the same sign. 

In the small x expansion at LO the quark anomalous dimensions vanish, hence the 
quark sector may be neglected, and the gluon evolution equation is found using the splitting 
function Pgg eq. ( p.lO|) in the Altarelli-Parisi equation eq. (| 



(3.21) 



n=0 •^—S.o J—£.o J—£.o 

where we have used the one loop form of as{t) as appropriate for a LO calculation, the 
coefficients are given by eq. (|3.15|) , and 7 is as in eq. ( |2.23| ). Since we retain only singular 
contributions to the splitting functions the lower limit ^0 = the integrations on 

the right hand side can be consistently set to zero. Differentiating l3oth sides w.r. to ^ this 
can again be cast in the form of a wave equation 



d^dC 

71=1 



^0 

00 A />.\r, 



In — 1 



0, 

(3.22) 
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which, when only the first singularity is retained, reduces to eq. (|2.22| ) with S = (the 
term contributing to S is NLO in the small x expansion). 

The quark does not evolve at LO, hence, in order to determine the evolution of F2, 
we must go to NLO. This can be done in complete analogy to the derivation of the double 
scaling evolution equations eq. ( |2.22|) in the previous section: in fact the derivation there 
can be viewed as a simplified version of a NLO calculation in the small x scheme, where 
only the (order ctg) contributions with n = 2 — m to the anomalous dimensions eq. ( |3.9| ) are 
retained. Thus, we first diagonalize the anomalous dimension matrix, which up to NLO 
has eigenvalues 



(3.23) 



The corresponding eigenvectors are given by 

Q+ = a^G+ + 0{a') Q. = -^G. + 0{a). (3.24) 
199 hi 

where the anomalous dimension is expanded as in eq. ( |3.12D . Notice that the large eigen- 
vector condition does not depend on the unknown NLO gluon anomalous dimensions, 
and the small eigenvector condition has still the simple form of eq. ( |2.24D , thanks to the 
color-charge relation eq. ( |3.16D . 

Transforming to x space, the large eigenvector component G+ is seen to satisfy 
eq. ( pi. 221 ), plus 0{as) corrections, which lead to the d term of eq. ( ^^.22] ) and to 0{as) 
corrections to the coefficients a„ in eq. ( |3.22| ). The quark equation is then found differ- 
entiating the large eigenvector condition wih respect to C and using the equation for 
with the result 

dC 9 L Jo {n-iy. J 

where the coefficients re given by 

a^ = {Ain4\n2CA/7rr''- (3-26) 

Although these evolution equations cannot be solved in closed form, a solution can be 
developed perturbatively in the usual way by noting, as we did at two loops (|2.41|) , that 
the solution (|2.27|) for G+ now acquires an extra set of terms 

~ J-J- Jo Jo Jo 

(3.27) 

on the right hand side. The full solution may then be developed iteratively by substitution 
of lower order solutions, leading to a power series expansion in terms of integrals over Bessel 
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functions, which can then be evaluated numerically. A similar expansion may be obtained 
in Mellin space by integrating explicitly the singlet renormalization group equation ( |2.7|) , 
and then performing the inverse transform (|2.12|) by choosing a contour which encircles 



the singularities in the region |A^| < 

For present purposes, however, we derive some simple analytic estimates which demon- 



strate analytically the nature of the solution [|2^ . The key observation which dictates the 
structure of the solution is the simple property of Bessel function z'^In-i{z) — (^z'^In(z)^ . 
This implies that if we determine iteratively the solution to the full gluon equation eq. ( |3.22| ) 
by then eq. ( |3.22| ) with the perturbation eq. (|3.27| ) takes at lowest order of the iteration 
the approximate form 

^ l'G{C, + ^ E an+1 (^^) ^n(27a). (3.28) 

But all Bessel functions in this expansion depend only on the scaling variable a and have 
the same asymptotic behaviour eq. (|2.26|) (which only sets in more slowly for higher values 



of n). Since p^In{'2'fa) is bounded above by $,^Io{2'ya) it follows that double scaling will 
always set in asymptotically provided the series Xl^i o^n+i (^^^^^)"^ converges uniformly, 
that is provided ^ < x^c)' '^^^^ much wider region than that in which double scaling 
would hold if we were to impose a "hard" boundary condition of the form at ^ = 0, 
namely ^ < x7o) ' because X{() rapidly falls as C increases. Such a boundary condition is 
precisely that which the simple argument based on the location of the singularity in 7^ 
eq. ( p.3|) would suggest. 

To see more clearly how such a power-like behaviour could arise, consider the evolution 
equation eq. ( |2.22| ) in the region of extremely small x, where the higher order terms give 
the dominant contribution to the series on the right hand side. Because the series ( |3.14| ) 
has unit radius of convergence it follows that lim """""^ — > 1. But then, setting a^+i ~ 

in the sum in eq. (|3.22| ), shifting the summation index by one unit, and using eq. ( |3.21|) , 
we have 

J2ar,+^X,{Cr ^6 • ■ • / d^' Gi^' X) = UO^- (3.29) 

Po ^ Jo Jo Jo dC 

Hence, asymptotically as ^ ^ 00 the evolution equation ( |3.22| ) becomes simply 

g^^GiiX) - A.(C)^ = I'GitC), (3.30) 

and the summation of all leading singularities leads to a damping term in the wave equation. 

We can now see how the behaviour obtains in the fixed-coupling limit: if the 
coupling is frozen, then A is just a constant, and the solution to eq. ( |3.3C1| ) is given in terms 
of the (double scaling) solution 6*0(^7 C) to the original wave equation 

G(e, = e^«Go(e, = x-^Goi^. 0- (3-31) 
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However there is no reason to fix the couphng. A solution with running couphng can be 
derived by the saddle point method , and turns out to give the double scaling behaviour 
(up to a small correction) 



2 



G(^, C) ~ ^e'^-Hm-mp- ^ (3.32) 



throughout the region ^ -C (a(o)_a)^ C^- When ^ is extremely large the power-like behaviour 



1^ i 



is found instead. This has essentially the "hard" form of eq. ( |3.31| ), with the large value 
of A evaluated at the initial scale, but it is confined to the extremely small region ^ ^ 

2 

A(o)2 ^ • -'-^ f^^^ ^^^^ ^^^^^ overestimate of the region where the power-like behaviour 
should arise: approximating the evolution equation with eq. ( p.30|) ignores the fact that 
the asymptotic behaviour of the coefficients only sets in rather slowly, and, by the time 
it does, an « ai. If this effect is taken into account the power-like region is further 
reduced. 

Even though of course all these estimates are only based on the asymptotic form of the 
coefficients a^, and thus they will not accurately reproduce subasymptotic and specifically 
subexponential corrections (such as the factor of ^ in eq. ( |3.32| )) and normalizations, they 
do correctly estimate the leading behaviours and their region of validity, as confirmed by a 
full numerical analysis. This allows us to conclude that, even though a power like behaviour 
of the form eq. ( |3.31j ) is generated very close to the C = boundary, it very rapidly dies off 
due to the running of the coupling, rather than spreading in the whole ^ > region as 
it would do if it were input to LO evolution. Furthermore, any hard power-like boundary 
condition will be hidden by this rise (unless it is even stronger than it): since A(0) is very 
large for reasonable choices of the starting scale (for instance A(4 GeV^) ~ 0.8) this means 
than scaling will appear very rapidly for all boundary conditions, except unreasonably 
hard ones such as with A ^ —0.8 at = 4 GeV^. 

Turning finally to the quark equation, we may take advantage of the fact that, through- 
out most of the (^, Q plane, (j(^, C) — A/'/o(27a): substituting this in the r.h.s. of eq. (|3.25|) , 
and neglecting the small eigenvector we get 

^ n=0 ^ ' 

2 2 

It follows that when ^ ^ ttCc^'', Q(^, C) still scales; for ^ ^ y^-C^^^, if we set a„ = 1 for all 
n, (|3.34| ) may be simplified yet further by using the relation Yl^oo^^-^n{z) = exp [^z{t + 
t~^)) , to give 

Qit ~ ArCe^^°+^'^/^°. (3.35) 



Q(^, ^ -f t'AT y~aj /.(27a). (3.34) 
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The quark anomalous dimensions may thus stiU produce a growth of F2 close to the 
boundary, which however again does not spread in the Q plane due to the running of the 
coupling. The actual size of the region where this growth will appear depends on the size 
of the coefficients a„. This, in turn, strongly depend on the choice of gluon normalization 
which, as discussed above, can be modified by changing factorization scheme. The latter 
could be fixed using momentum conservation if we knew the NLO small x gluon 
anomalous dimensions. Since we do not, the effect of the quark anomalous dimensions on 
the relative normalization of quark and gluon distributions (and thus on the relative size 
of F2 and Fl) turns out to be strongly dependent on the choice of factorization scheme: 
while the corresponding to the MS or DIS schemes [^| give large (but very different) 
effects, the Qo-schemes [^, being less singular, lead to substantially smaller effects, while 
in some schemes j^^j there is almost no effect at all except at very small x very close 
to the boundary. 

These analytic results are all closely supported by detailed numerical investigations, 
either retaining only the singular terms in the anomalous dimensions [p6|,|32[|, or including 



also the full one and two loop contributions in the double leading scheme |jT^,0 or using 
some other procedure p5|] . When care is taken to consistently factorize residual nonper- 
turbative effects into the boundary conditions, all these analyses reach essentially the same 
conclusions. The leading (Lipatov) singularities have only a negligible effect throughout 
the measured region, so small that there is still no empirical evidence for them, not only 
because they don't affect the shape of F2 in most of the (^, Q plane, but also because the 
coefficients are so small. The subleading quark anomalous dimensions have little effect 
of the shape of F2, explaining the success of conventional perturbation theory and double 
scaling p6| , p!9[1 . However they can have a substantial effect on the scale Qo at which soft 

and on 



initial distributions are input, raising it to around 2GeV in the MS scheme p6, 25 



the relative normalization of quark and gluon, and thus on the size of F^ as deduced from 

However these latter effects are strongly factorization scheme 
and this scheme dependence could only be 
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the measured F2 

dependent, both in sign and magnitude 
resolved theoretically by a complete calculation of the subleading singularities of the gluon 
anomalous dimension, or phenomenologically by a direct measurement of F^. 



4. Outlook 

The observation of double asymptotic scaling at HERA is a striking success of a 
perturbative QCD prediction made now more than 20 years ago 0. The explanation of 
the effect turns out to be somewhat subtler than it might seem at first: even though it is a 
direct consequence of the singularity structure of the leading-order perturbative evolution 
equations, its stability in a region where higher order effects might naively be expected to 
be important is due to the all order cancellation of double logarithmic singularities, the 
unexpected smallness of the coefficients of the remaining single logarithmic singularities, 
and then finally the reduction in their impact when the effect of the running coupling is 
included in a way consistent with the renormalization group. 

This result has interesting ramifications from both the phenomenological and the 
theoretical point of view, which have just started to be explored. The success of NLO 
perturbation theory in the HERA region strongly suggests that this may be an ideal place 



25 



to perform precision tests of perturbative QCD, since most of the poorly known low-energy 
effects (such as higher twists) are either absent or negligible here. For example, can 



be measured from existing HERA data [|T9[ with a precision comparable to that of all 



other existing deep-inelastic experiments combined. As more data in the very small x, 
not-so-small region become available, they may shed light on subtle issues of scheme 
dependence, and provide information on the behaviour of parton distribution which are 
input to perturbative evolution. On the more theoretical side, it would be desirable to 
derive perturbation theory at small x in a way which treats the two large scales symmetri- 
cally from the outset, rather than solving renormalization group equations with respect to 
one scale while including the summation of the other scale in the anomalous dimensions, 
as we did here. This, besides being interesting for its own sake, could shed light on the 
dynamics of perturbative QCD in the high energy regime. 
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Figure Captions 



Fig. 1. Scaling plot of ttie experimental data on the proton structure function F2(x, Q^) 
from (a) the 1992 HERA run (adapted from ref. 0) , (b) the preliminary 
analysis of the 1993 HERA run [| (adapted from ref. [|), (c) the 1993 run 0], 
(d) the preliminary analysis of the 1994 run (HI only) (adapted from ref. [^). 
The diamonds are ZEUS data and the squares are HI data. Only points with 
p, a > 1 are included in the plots b)-d) , and in (d) only those with > 5 GeV^ 
(see text). The straight line shown is the predicted asymptotic double scaling 
behaviour (with fitted normalization). 

Fig. 2. The (^,C) plane, showing the backward light cone at the point {$,',C): curves of 
constant a (the hyperbolae) and lines of constant p. The origin is chosen in 
such a way that the small x approximation of the evolution equations is valid for 
positive ^, and perturbation theory breaks down for negative ( where ctg grows 
too large. The hatched area x ^ exp(— as(to)^/«s(^)^) indicates the region 
where parton recombination effects are expected [|11| to lead to breakdown of 
simple perturbative evolution due to higher twist corrections. 

Fig. 3. Double scaling plots of i^F-^l' against i) a and ii) p, (a) with rescaling of the LO 
asymptotic behaviour eq. (|2.36|) ; (b) with NLO rescaling eq. ( |2.44|) . The stars are 



preliminary low (l.SGeV^ < < 15GeV^ 3.5 10"^ < x < 4.0 lO'^) points 
from ZEUS 

Fig. 4. The terms summed in the various expansions of the anomalous dimensions and 
associated splitting functions: a) the standard (large-x) expansion ( |3.6|) , b) the 
small-x expansion (|3.9| ), and c) the 'double- leading' expansion ( |3.11|) . Leading, 
sub-leading and sub-sub-leading terms are indicated by the solid, dashed and 
dotted lines respectively; m denotes the order in a, while n denotes the order in 
1/A^. Singular terms are marked as crosses, while terms whose coefficients known 
at present (for 7^) are marked by circles: the term which leads to double scaling 
is marked with a star. 
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